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We study the polarized Bjorken sum rule at low momentum transfer squared < 3 
GeV^ in the twist-two approximation and to the next-to-next-to-leading order accuracy. 
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'p^; 1 Introduction 

The spin structure of a nucleon is one of the most interesting problems to be resolved 
' within the framework of (nonperturbative) Quantum Chromodynamics (QCD). In par- 

ticular, the singlet part of the parton distribution functions 

> : / 

00 

' where / is a number of active quarks, is intensively studied, because there is strong 

. disagreement between the experimental data for its first Mellin moment and corresponding 

theoretical predictions. This disagreement is usually called a spin crisis (see, for example, 
reviews in [1]). 

Here we consider only the non-singlet (NS) part, which the fundamental Bjorken sum 
rule (BSR) holds for [2] 



dx. 



rr"(Q2)= / k(x,Q2)-g?(x,Q^ 
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It deals with the first moment (n = 1) of NS part of the structure function (SF) gi{x, Q"^)- 
For the case n = 1, the corresponding anomalous dimension of Wilson operators is zero 
and all the Q^-dependence of T\~^{Q'^) is encoded in the coefficient function. 
Usually, BSR is represented in the form 

CO P~'^ ( 

rr'^(Q^) = f^^5(Q^) + E^%iS^> 

where the first term in the r.h.s. is a twist-two part and the second one is a contribution 
of higher twists (HTs). 

At high values the experiment data for rj'~"(Q^) and the theoretical predictions [1] 
are well compatible with each other. Here we will focus on low values, at which there 
presently exist the very precise CLAS [3, 4] and SLAC [5] experimental data for BSR. 
On the other hand, there also is a great progress in theoretical calculations: recently, the 
terms ~ are evaluated in [6]. 
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Figure 1: (color online). CLAS [3, 4] and SLAC [5] experimental data for BSR and < 3 
GeV^. The curves represent theoretical predictions obtained in the analyzes carried out by two 
groups: Khandramai, Pasechnik, Shirkov, Solovtsova, and Teryaev (KPSST) [7] and Kotikov 
and Shaikhatdenov (KS). 



2 Basic formulae 

In our analysis we will mostly follow the analyses done by the Dubna-Gomel group [7, 8]. 
We try, however, to resum the twist-two part with the purpose of reducing a contribution 
coming from the HT terms. 
Indeed, there is an interplay 

• between HTs and higher orders of perturbative QCD corrections (see, for example, 
[9], where the SF xF^ was analyzed). 

• between HTs and resummations in the twist-two part (see, for example, application 
of the Grunberg approach [10] in [11] to the study of SFs F2 and Fl) 

The twist-two part of BSR has the following form (see, for example, [7]) 

i^V5(Q') = l-4A(Q2), (1) 

where the term A{Q'^) looks like 

A(Q2) = a.(Q2)(i + f;c7,aJ(Q2)) (^„^(q2) ^ . (2) 

The first three coefficients Ci, C2 and C3 are already known (see [6, 12] and references 
therein) . 

We will replace the above representation (2) by the following one 



where 

oo 

A(Q2) = a,(Q2)(i + ^ Cua};{Q^)) (4) 
fe=i 

and Ck can be obtained from the known C^: 

C'i = Ci+4, 6-2 = C2 + 8Ci + 16, (72 = C3 + 8C2+4C! +48Ci + 64. (5) 

The reason behind this transformation is as fohows: the CLAS experimental data [3, 4] 
demonstrate that V\~^{Q'^ — > 0) — >■ 0. Therefore, in the case when the HT corrections 
produce small contributions at — >■ ^ we see that 

Ens{Q^ ^^)^^. (6) 

Since the strong coupling constant a.^lQ'^ — > A^) — > 00, it is seen that the form (3) 
behaves much like the CLAS experimenatal data. Indeed, 

Ens(Q^^^^) = ^ ^0. (7) 

^ 1 + 4A(Q2 ^ A2) ^ ^ 

As Aq(-,£, ~ 0.01 is rather small, one can conclude that the above representation (7) 
agrees with experiment at very low values. 

Note, however, that we have a very small coefHcients of A{Q^) and A{Q'^). Thus, for 
small but nonzero values the above representations (1) and (2) lead to similar results 
(see Fig. 1, where we restricted our consideration to the next-to-next-to-leading order 
(NNLO) accuracy). As is seen in Fig. 1, the theoretical predictions are not too close to 
the shape of the experimental data. 



3 Grunberg approach 

At ~ 0, the value of the strong coupling constant is very large. Thus, in our approach 
it is better to avoid the usage of series like 

00 

J2C\a';iQ'), (8) 

k=l 

as in Eqs. (2) and (4). 

Instead, it is convenient to use the Grunberg method of effective charges [10], i.e. to 
consider the variables A((5^) and A((5^) as new effective "coupling constants" ^, which 
have the following properties: 

• shifted arguments 

for the variables A(Q^) and A.{Q^), respectively, with 

i^fc = e^i/^°, Dk = e^'/^°, (10) 

which are in turn responsible for the vanishing of the coefficients Ci and Ci in a 
series similar to (8). Moreover, these shifted arguments (9) provide also a strong 
reduction in the magnitudes of the coefficients Ck and {k >2). 

• new Pi {i > 2) coefficients of the corresponding ^-functions, which are responsible 
for the vanishing of the coefHcients Cfc and Ck (fc > 2). 

^It is for sure qucstionalolc, but in the KPSST analysis [7] /i4 ~ at the next-to-next-to-next-to-leading 
order (N-^LO) accuracy. 

similar application can be found in the recent paper [13]. 
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Figure 2: (color online). The experimental data are as in Fig.l. The curves show theoretical 
predictions based on Eqs. (1) and (3), which is called a KPSST-like analysis [7] and a KS 
one, respectively. For all cases the Grunberg approach [10] is used with a standard coupling 
constant and with a "frozen" one, when a = 1.5. 

However, a straightforward application of the Grunberg approach to the variables 
A((5^) and A(Q^) is not as convenient, because the coefficients Ci and Ci are positive 
and the values are very small. It is in contrast with its direct applications, where the 
coefficients Ci and Ci are negative [14] and/or the values are not so small [11, 15]. 

So, the new arguments Q^/D^ and Q^/D^ have now very small values and, as a result, 
we have to use the Grunberg approach associated with something else. One of the ways 
is to use a so-called "frozen" coupling constant. 

4 "Frozen" coupling constant 

We introduce freezing of the coupling constant by altering its argument = 
+ aMp, where Mp is a p-meson mass and a is some free parameter (usually, a = 1 was 

used. See, for example, [16]). 

Thus, in the formulae of the previous sections the following replacement should be 

done (a list of references can be found in [17]): 

asiQ^)^afr{Q^) = asiQ^ + aM^p) (11) 

In the analyzes given below we restrict ourselves to the next-to-leading order (NLO) 
(i.e. a^) approximation. The consideration of two even higher order corrections is in 
progress. 

The cases with a = 1.5 and a = 2 are shown in Figs. 2 and 3, respectively. It is seen 
that the best agreement with experimental data is achieved in the case of representation 













r I 








KPSST(as^) 




KS(as^) 






KPSST(frozen ttg^) 




. 1 . . , , 1 , 


KS(frozen a-^) 



Modified via A 



0.18 



0.12 - 



0.06 












KPSST(as") 
KS(as^) 
KPSST(frozen a^^) 
KS(frozen a^^) 



0.5 



1.5 

(GeV^) 



2.5 



Figure 3: (color online). As in Fig 2 but a = 2. 



(3) and a = 1.5. Also, the standard form (1) and theoretical predictions obtained with 
a = 2 are well consistent with each other. 

In contrast with the analyzes carried out with a standard coupling constant (see Fig. 
1), we observe that the shape of theoretical predictions and the form of experimental data 
are close enough to each other at small values. 



5 Conclusion 

The analysis of the Bjorken sum rule performed within the framework of perturbative 
QCD is presented at low Q^. It features the following important steps: 

• The new form (3) for the twist-two part was used. It is compatible with the observa- 
tion Ens{Q'^ ^ 0) ^ 0, coming from the experimental data (if HTs are negligible). 

• The application of the Grunberg method of effective charges [10] in a combination 
with a "frozen" coupling constant provides good agreement with experimental data, 
though with a slightly larger freezing parameter (1.5M^ instead of M^). ^ 

Further elaborations to be undertaken include taking into account the and cor- 
rections to our analysis, as well as the study of HT corrections and their correlations with 
a freezing parameter a (in front of M^). We also plan to add to our analysis an analytic 
coupling constant [18], which has no the Landau pole and leads usually to the results, 
which are similar to those obtained in the case of the "frozen" coupling constant [17, 19]. 

A.V.K. thanks the Organizing Committee of 20th International Symposium on Spin 
Physics (SPIN2012) for invitation and support. 



■^It seems that the value of the parameter a depends on the order of perturbation theory. We plan to study 
this dependence in our forthcoming investigations. 
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